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A combination of classical plate theory and a supersonic aerodynamic model is used to 
study the aeroelastic flutter behavior of a proposed thermal protection system (TPS) for the 
NASA HIAD. The analysis pertains to the rectangular configurations currently being tested 
in a NASA wind-tunnel facility, and may explain why oscillations of the articles could be 
observed. An analysis using a linear flat plate model indicated that flutter was possible well 
within the supersonic flow regime of the wind tunnel tests. A more complex nonlinear 
analysis of the TPS, taking into account any material curvature present due to the restraint 
system or substructure, indicated that significantly greater aerodynamic forcing is required 
for the onset of flutter. Chaotic and periodic limit cycle oscillations (LCOs) of the TPS are 
possible depending on how the curvature is imposed. When the pressure from the base 
substructure on the bottom of the TPS is used as the source of curvature, the flutter 
boundary increases rapidly and chaotic behavior is eliminated. 
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plate length 

modal coordinate of Kevlar/general plate 
speed of sound 
plate width 

modal coordinate of Nextel 

Fourier coefficient for the imposed-pressure geometry 

damping matrix 
plate bending rigidity 
Young’s modulus 
Airy stress function 

modal coordinate of the Airy function expansion 

rise height of a curved plate 
plate thickness 
spring constant 
stiffness matrix 
mass matrix 
mass/area of plate 
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= static in-plane loading (tension) in the x-direction 
= static in-plane loading (tension) in y-direction 
= static in-plane loading (tension) in the x-y-direction 

= eigenvalue for the flat plate flutter problem/imposed pressure 
= aerodynamic pressure 
= eigenvector for the flat plate flutter problem 
= generalized force 

= free stream velocity 

= plate deflection 
= curved plate rise function 
= mode shape for the Airy function 

= mode shape of the Kevlar 
= mode shape of the Nextel 
= flow density 

= non-dimensional plate length 
= non-dimensional plate width 
= structural damping ratio 
= Poisson’s ratio 

= fundamental natural angular frequency of a plate 


I. Introduction 

T he purpose of the following analysis is to characterize the flutter behavior of the proposed thermal protection 
layer (TPS) on the NASA Hypersonic Inflatable Aerodynamic Decelerator (HI AD). The TPS is composed of 
several layers of different materials, each with unique physical properties. At certain flight conditions, characterized 
by flow velocity and density, interactions between aerodynamic forces and the structure may cause a dynamic 
instability known as aeroelastic flutter. Since this instability could damage the TPS and possibly expose the HIAD to 
the high temperature flow of re-entry, it is important to determine the conditions in which flutter may occur. 

While the aeroelastic behavior of the HIAD of great interest, its complex geometry makes it difficult to analyze, 
especially with respect to flutter. Thus, as a first step, the present analysis is based on the current NASA wind-tunnel 
tests of rectangular TPS articles 1 , which are mounted on metallic tiles and exposed to hypersonic flow. A flutter 
analysis of this geometric shape is less complex, and it may help explain why oscillations could be observed during 
testing. While flutter is a possibility, excitation of a base structural component or turbulent boundary layer 
fluctuations cannot yet be ruled out as sources of oscillation. In any case, this investigation is an important 
foundation for the eventual consideration of the aeroelastic behavior of the TPS in its full flight configuration. 

II. Theoretical Models 

Several structural models have been developed for the rectangular test article, including both flat and curved 
plate models. The test article is a multi-layer structure consisting of Nextel 400-BF20, Pyrogel 6650, and 
Aluminized Kapton Kevlar Laminate. Curvature was estimated based on photographs of the test samples in an 
earlier AIAA report by Hughes, et al 2 . In all cases, piston theory is used for the aerodynamic pressure. This linear 
unsteady theory assumes the induced aerodynamic pressure is a result of the instantaneous plate deflection only, and 
thus memory effects are ignored. 

A. Linear Flat Plate Model 

The first proposed structural model is a two-dimensional linear system consisting of two flat plates separated by 
a spring-damper, shown in Figure 1 . The Nextel and Kevlar are the flexible plate members while the inner Pyrogel 
layer is a spring-damper. Since the Pyrogel is relatively thick and sponge-like, it can be compressed as a result of 
deflection of both the Nextel and Kevlar. The model assumes that the Nextel and Kevlar are connected to the 
Pyrogel at every point along the length of the sample, which is an approximation whose validity depends on how the 
layers are stitched together. 
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In-plane loading (tension) will also be considered, since some tension may be applied in the wind-tunnel test. A 
uniform uni-directional supersonic flow along the plate length with velocity U is also present. Simply supported 
boundaries are shown here; however, other boundary conditions such as clamped and free can be imposed and will 
be discussed later. Note that by setting the spring stiffness and damping ratio of the Pyrogel to zero, the system 
decouples and reduces to a single Nextel (or Kevlar) plate which may also be of interest. By increasing the bending 
stiffness of the Kevlar to a large value, the system becomes a single Nextel plate on a Pyrogel foundation. Flutter 
points and natural frequencies were calculated for all of these configurations and are discussed in the results section. 

U 

► 


Nextel 

Pyrogel 

Kevlar 



Figure 1. Flat plate model of the TPS system with simply supported 
boundary conditions and flow velocity U. 


1. Governing Equations 

A classical Lagrangian energy approach was used to derive the equations of motion for this system. The general 
formulation is provided below, and the detailed expressions for the modal equations can be found in the appendix. 
The total potential energy of the system is: 
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( 1 . 2 ) 

Note that only bending potential energy is used here. An assumption has been made that the plate layers do not 
sufficiently stretch during deflection, which is the case for flat plates with oscillations substantially less than the 
thickness 3 . For larger deflections, in the case of initial/imposed curvature or post-flutter limit cycle oscillations, the 
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plate material may stretch, which would require the addition of a stretching potential energy. This addition will be 
important for the curved plate models in the next sections. The total kinetic energy for this system is: 
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The deflection of both the Kevlar and Nextel can be expanded using the following modal series: 


(1.3) 
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The modal functions (p n and y/ r are the products of beam modes that depend on the boundary conditions of the 
problem. Only one mode is retained in the spanwise direction, which is a valid assumption since the plate is flat and 
the flow is only along the plate length. Substituting Eq. (1.6) into Eq. (1.5) results in a Lagrangian in terms of the 
two modal coordinates a n and b r . Lagrange’s equations for the problem are: 
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The generalized force terms Q n and Q r incorporate the aerodynamic pressure and any structural damping that 
may be present. In this analysis, the aerodynamic pressure is given by the piston theory approximation for 
supersonic/hypersonic flow: 


A P = PA 


dw + dw 
00 8x dt 


(l.B) 


The aerodynamic pressure is only applied directly to the outer Nextel, so the Kevlar is not exposed to the flow. 
The stability of the system can be studied by further expanding the time dependent modal coordinates as follows: 
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(1.9) 

b r =b r e pt 

The variable p is a complex number which characterizes both system damping and frequency. Substituting Eq. 
(1.9) into the equations of motion derived from Eq. (1.7) produces the following eigenvalue problem for stability of 
the system: 


(p 2 M+pC + K){q} = { 0} 


( 1 - 10 ) 


The quantities M, C, and K are the mass, damping, and stiffness matrices, respectively. This eigenvalue problem 
is solved using the P-K method, whereby the eigenvalues are calculated for a given frequency until convergence is 
attained, and then stability is determined by observing any sign changes that occur in the damping, all the while 
increasing a flow parameter in small increments. 


B. Parabolic Curvature Model 

This three-dimensional nonlinear structural model accounts for any small amount of curvature that the TPS may 
exhibit when tensioned across the rigid tile. The test configuration from which this model is based is shown in 
Figure 2, and the model representation is given in Figure 3. Visual inspection of the test article indicates some 
(approximately) parabolic curvature due to the edge restraint bracket. The model assumes there is a space between 
the center of the bottom Kevlar layer and the center of the tile which is equal to the curvature rise height (H). The 
curvature is assumed to be a parabolic function in both the streamwise and spanwise directions, although results will 
also be computed for the case with streamwise curvature only. The maximum rise height of the tensioned TPS will 
be estimated as on the order of five times the thickness of the Nextel, which is approximately 0.5 cm. Calculations 
will also be performed for a range of smaller rise heights, since interesting flutter behavior may arise in this regime. 
Note that the flutter amplitudes of curved plates are on the order of the rise height, rather than the plate thickness 4 , as 
is the case for flat plates. This is why nonlinear large -deflection theory is required. 

Additionally, the piston theory aerodynamic pressure in Eq. (1.8) must be modified to account for the presence 
of curvature, resulting in an additional static loading term proportional to the initial plate slope. This static load is 
critical in determining the flutter boundary. When the curvature is set to zero, this model reduces down to a three- 
dimensional nonlinear version of the previous two-dimensional linear flat plate model with standard piston theory 
aerodynamics. 



Figure 2. TPS test article tensioned across rigid metallic tile, 
exhibiting streamwise and spanwise curvature 2 . 
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Figure 3. Theoretical structural model of the TPS article with initial 
parabolic curvature. 


1. Governing Equations 

The equations of motion for this system were derived using the Airy stress function with stress boundary 
conditions for the in-plane deformation. The bracket shown in Figure 2 may apply stress at the boundaries, but it is 
likely smaller than the stress due large amplitude out-of-plane bending. The bracket also cannot enforce a condition 
of vanishing slope at the boundaries, since a small space exists between the material and the bracket. These two 
observations imply stress-free boundary conditions, although the effect of applied stresses will be considered as 
well. The out of plane boundary conditions will be approximated as simply-supported, since some degree of rotation 
is possible within the space between the bracket and TPS. This configuration will be a lower-bound for flutter, since 
simply supported edges give rise to a lower flutter boundary than clamped edges (see Table 7). 

The governing equations for the TPS with parabolic curvature were derived using the same Lagrangian approach 
for the flat plate analysis, with some modifications to each individual plate layer. The total potential energy for a 
single plate is now: 


u l - plate =u B +u s 

The bending potential energy is the same as Eq. (1.2), and the stretching potential energy is: 


( 1 . 11 ) 
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Where F is the Airy stress functions defined by: 
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The quantities N x , N y , and N xy are the in-plane stress resultants. The Airy function also satisfies the following 
equation of strain compatibility, which is modified for the case of initial curvature: 
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The function z is the parabolic plate rise function defined by: 
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The modification to the piston theory aerodynamic model, to account for plate rise, is given by: 
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Eq. (1.14) along with Eq. (1.12) and Eq. (1.11) comprise the modified energy formulation of Von Karman’s 
large deflection theory. The solution procedure begins by again expanding all time-dependent functions into modal 
series: 


w(x, y,t) = Y j a n (t)</> n (x, y) 
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Again, both modal functions must satisfy the boundary conditions. The particular solution F is the solution to 


the homogenized version of Eq. (1.14), and is associated with the static loading of the plate : 

F _ , Nj 


-Nxyxy 


(1.18) 


The quantities N x , N y , and N xy also satisfy stress boundary conditions. For stress free boundaries, these 
quantities are zero. When stress is imposed at the boundaries, it is assumed that no shear stress is applied. A solution 
to the full strain compatibility condition is obtained via Galerkin integration, which allows the stress coordinate /. 

to be written in terms of the deflection coordinate a n . The detailed equations of motion for a single plate are given 
in the appendix. 


C. Imposed Pressure Model 

In this analysis, the shape of the TPS is assumed to be the same as the metallic tile substructure. Unlike the 
previous model, this configuration has a flat middle surface with parabolic curvature only at the edges. The 
experimental configuration from which the model is based is shown in Figure 4, and the theoretical model is shown 
in Figure 5. This model also assumes that when the TPS is secured across the metallic substructure, a static pressure 
is imposed on the Kevlar. This is the pressure required to generate the initial geometry of the TPS, which can be 
solved via a static inverse problem. This model can also be interpreted as a three-dimensional, nonlinear version of 
the flat plate model with an additional generalized force due to the imposed pressure. Note that the infinite rigidity 
of the tile substructure will not be considered. That would require enforcing the condition that there be no deflection 
at every point where the TPS contacts the tile. 
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Figure 4. TPS test article tensioned across rigid tile, exhibiting a 
flat middle surface and parabolic edge curvature 1 . 


U 



Pyrogel 


Figure 5. Curved TPS model due to imposed pressure from the 
rigid tile. 


1. Governing Equations 

A similar nonlinear stress formulation from the previous curvature model is used here, but with several key 
modifications. Again, these modifications apply to a single plate, but the equations for full multi-layer TPS can be 
easily derived following the procedure for the flat plate model in section A. In this case, the rise height of the initial 
curvature function is set to zero, such that the condition of a flat plate is attained. The curvature is imposed via a 
generalized force, which is calculated in the following manner. The desired curvature functions (in the x direction 
only) are: 
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Where jq and x 2 are the positions along the plate length where the edge curvature begins and ends, respectively. 
Since measurements of the rigid base are not yet available to the authors, the region of curvature at the edges was 
estimated to be ten percent of the length of the plate (based on the presentation of the curvature in Figure 4). For the 
purposes of convenient numerical integration, Eq. (1.19) is expanded in a Fourier series of the following form 

Z=Y, C n a n^y) C 1 ' 20 ) 


The Fourier coefficients c n are given in the appendix. Not that while only a single index is written, summation is 
in both x and y directions. The resulting surface from the Fourier series expansion is shown in Figure 6. 



Figure 6. Fourier series representation of the rigid base geometry. Panel length and width 
have been non-dimensionalized. 


The imposed pressure p is calculated using the static form of Von Karman’s large deflection equations: 
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Numerical integration of (1.21) results in the generalized force term that can then be used in the dynamic 
equations of motion. The full expression is provided in the appendix. Again, stress-free and imposed-stress in-plane 
boundary conditions will be considered with simply supported out of plane conditions. 


III. Material Properties 

The physical properties for the three TPS materials are given in Table 1, while the nominal geometry and 
mechanical properties are given in Table 2. The materials studied in this paper include Nextel 400-BF20, Pyrogel 
6650, and Aluminized Kapton Kevlar Laminate. 
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Table 1. TPS material properties from material documentation and various literature sources 5 ' 7 


Material 

2x Pyrogel 6650 

lx Kevlar 

2x Nextel 440* 

Thickness - h (m) 

0.00598 

0.000127 

0.001016 

Density - p (kg/m 3 ) 

112 

1440 

3050 

Mass/area - m (kg/m 2 ) 

0.67 

0.18 

3.1 

Young’s (tensile) modulus - E (Pa) 

3.80xl0 6 

1.13xl0 n 

1.90xl0 n 

Poisson ratio - v 

0.20 

0.36 

0.36 

Bending rigidity - D (Pa m 3 ) 

- 

0.022 

2.39 

Spring stiffness - k (Pa/m) 

6.60x10 s 

- 

- 


Table 2. Nominal configuration geometric and structural properties 


Sample Length - a (m) 

0.1778 

Sample Width - b (m) 

0.127 

Nextel Damping Ratio - 

0 

Kevlar Damping Ratio - 

0 

In-plane loads (tension) - N x , N v (N/m) 

0 


*Note: The bending rigidity for the Nextel was calculated assuming two sandwiched layers as in the layup 
configuration 1 ; however, it was determined from preliminary experiments that the bending rigidity of two layers was 
not equal to twice that of one layer. In fact, the measurements for both cases were nearly the same. As a result, and 
for the purposes of determining the lower bound on flutter, the bending rigidity of one layer of Nextel was used in 
all calculations. Also, the two layers of Pyrogel are modeled as a single layer with twice the nominal thickness. 
Additional structural tests to measure the bending properties directly are currently in progress. 


IV. Natural Frequency Results and Analysis 

The natural frequency results are presented in the following section for all three theoretical models. For the 
linear flat plate model, calculations were made using an eigenvalue solver, while the calculations for the two 
nonlinear models were made using a time-marching algorithm in MATLAB. 

For each of the three structural models, various configurations of the TPS materials were investigated. The 
following specifications are used to reference to the different configurations: when a reference is made to “Nextel 
Alone” or “Kevlar Alone,” the other TPS layers are not included in the calculation but all other properties of the 
structural model remain the same. This would apply if there is no knitting connecting any of the layers together in 
the test setup. When a reference is made to a material on a “Pyrogel foundation,” the material rests on the Pyrogel 
layer which is itself resting on a rigid base. When all three layers of the TPS are included in the analysis, it will be 
referred to as the “3-layer TPS.” These specifications will be used in the discussion of the flutter results as well. 

A. Linear Flat Plate Model 

The natural frequencies were calculated for several configurations of the linear flat plate model, including the 
Nextel alone, Kevlar alone, and 3-layer TPS, with and without in-plane loads. In addition, both simply- supported 
and clamped boundary conditions were considered. All of these results are presented in Table 3-5. 


10 

American Institute of Aeronautics and Astronautics 





Table 3. Natural Frequencies of various Nextel configurations in rad/s, 


Mode 

Nextel 
Alone, 
Nominal 
(BC: SS) 

Nextel Alone 
N x =N y = 
1000 N/m 
(BC: SS) 

Nextel on 
Pyrogel 
Foundation 
(BC: SS) 

Nextel Alone, 
Nominal 
(BC: Clamped) 

Nextel Alone 
N x =N y = 
1000 N/m 
(BC: Clamped) 

Nextel on 
Pyrogel 
Foundation 
(BC: Clamped) 

1 

1405 

1694 

25353 

2641 

2845 

25451 

2 

2828 

3132 

25471 

4287 

4531 

25674 

3 

5201 

5511 

25842 

7019 

7284 

26269 

4 

8523 

8837 

26710 

10767 

11043 

27508 

5 

12795 

13110 

28363 

15491 

15774 

29677 

6 

18015 

18331 

31069 

21175 

21463 

33002 

7 

24184 

24502 

35009 

27849 

28141 

37635 

8 

31303 

31621 

40257 

35460 

35755 

43569 

9 

39371 

39689 

46806 

41580 

41859 

48679 

10 

48388 

48706 

54609 

48227 

48580 

54467 


Table 4. Natural Frequencies of various Kevlar configurations in rad/s 


Mode 

Kevlar 

Alone, 

Nominal 

(BC: SS) 

Kevlar 
Alone 
N X =Ny= 
1000 N/m 
(BC: SS) 

Kevlar on 
Pyrogel 
Foundation 

(BC: SS) 

Kevlar Alone, 
Nominal 
(BC: Clamped) 

Kevlar Alone 
N X =Ny= 
1000 N/m 
(BC: Clamped) 

Kevlar on 
Pyrogel 
Foundation 
(BC: Clamped) 

1 

323 

2289 

60554 

607 

2602 

60556 

2 

651 

3280 

60557 

986 

3647 

60561 

3 

1196 

4521 

60565 

1614 

4922 

60575 

4 

1960 

5916 

60585 

2476 

6366 

60604 

5 

2943 

7444 

60624 

3563 

7958 

60658 

6 

4143 

9111 

60695 

4870 

9698 

60749 

7 

5562 

10924 

60808 

6405 

11599 

60891 

8 

7199 

12893 

60979 

8155 

13663 

61100 

9 

9055 

15029 

61226 

9563 

14982 

61303 

10 

11129 

17340 

61567 

11092 

17846 

61560 


Table 5. Natural frequencies of the 3-Layer TPS in rad/s. 


Mode 

3-Layer TPS, Nominal 
(BC: SS) 

3-Layer TPS, 
N X =Ny= 
1000 N/m 
(BC: SS) 

3-Layer TPS, 
Nominal 
(BC: Clamped) 

3-Layer TPS, 
N X =Ny= 
1000 N/m 
(BC: Clamped) 

1 

726 

1001 

1364 

1566 

2 

1461 

1758 

2211 

2456 

3 

2680 

2992 

3608 

3877 

4 

4370 

4691 

5493 

5780 

5 

6492 

6824 

7792 

8097 

6 

50781 

50799 

50812 

50835 

7 

50818 

50856 

50882 

50927 

8 

50936 

51005 

51074 

51152 

9 

51219 

51331 

51488 

51611 

10 

51787 

51952 

52264 

52441 
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The natural frequencies for the various linear flat plate configurations have several important properties. First, 
the effect of applying tension at the boundaries of 1000 N/m had surprisingly little effect on the Nextel alone, while 
the natural frequencies of the Kevlar increased significantly. This same loading had a minimal effect on the 3 -layer 
system, which suggests that the stiffer Nextel governs the TPS dynamics and should also govern the flutter behavior. 
However, the natural frequencies of the 3 -layer TPS are quite different than the natural frequencies of the Nextel 
alone. The first five modes are lower than the corresponding modes of the Nextel, which is due to the added inertia 
of the Pyrogel and Kevlar. The next five modes are significantly higher than those of the Nextel alone, but this is 
because five modes were used for both the Nextel and Kevlar in the 3 -layer TPS eigenvalue calculation. Since the 
Nextel is orders of magnitude more rigid than the Kevlar, the Kevlar behaves as if there is a Pyrogel foundation, 
resulting in the higher natural frequencies. Finally, the clamped boundary conditions tend to increase the natural 
frequencies by approximately a factor of 2 with respect to simply supported boundaries for the lower modes. This 
effect diminishes when the Pyrogel foundation is added to the Nextel and Kevlar because the Pyrogel spring 
stiffness is quite large compared to the plate bending stiffness. 

B. Parabolic Curvature Model 

Although this model includes the Nextel, Pyrogel, and Kevlar, natural frequencies of the Nextel alone were also 
computed. Figure 7 illustrates the effect of streamwise and spanwise curvature on the natural frequencies of the 
Nextel alone. When H/h is zero, there is no curvature, and a flat plate is attained. The natural frequencies for this 
case are nearly the same as those calculated in the linear flat plate analysis. Two different computational/numerical 
methods were used for these calculations, so obtaining nearly the same results provides some degree of verification 
to both methods. 

When the Nextel has curvature in both directions, all of the natural frequencies increase as H/h increases (H and 
h are the rise height and plate thickness, respectively). With streamwise curvature only, shown in Figure 8, the 
natural frequencies intersect at various locations corresponding to different values of H/h. At these intersection 
points, one might expect the plate to be more susceptible to flutter since eigenvalue coalescence has occurred. 4 It 
will be shown later that in this case, the flutter boundary is usually lower than for a flat plate. 



Figure 7. First four natural frequencies of the Nextel alone as 
a function of H/h, with both streamwise and spanwise 
curvature. 
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H/h 

Figure 8. First four natural frequencies of the Nextel alone as 
a function of H/h, with streamwise curvature only. 



Figure 10. First four natural frequencies of the 3-layer TPS as 
a function of H/h, with both streamwise and spanwise 
curvature. 



Figure 9. First four natural frequencies of the 3-layer TPS as 
a function of H/h, with streamwise curvature only. 
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The natural frequencies for the 3 -layer TPS are given in Figure 9 and Figure 10. The general trend of these 
natural frequencies is similar to those for the Nextel alone, while the magnitudes of the frequencies are lower. In the 
case of streamwise curvature only, the points of intersection between the modes are nearly identical to the Nextel. 
Again, this would indicate that the flutter behavior of the Nextel is similar to that of the 3 -layer TPS, and the results 
in the next section indicate that this is generally the case. 

C. Imposed Pressure Model 

The natural frequencies for this model are the same as the parabolic curvature model with H/h equal to zero, 
corresponding to a flat plate. This is because the curvature is induced by an external forcing, rather than the 
geometric configuration of the material. 


V. Flutter Results and Analysis 

The flutter results that follow will be presented in terms of dimensional parameters for easier comparison to the 
supersonic/hypersonic wind tunnel tests. The Mach number will be fixed to a value of 5, and the flow density will 
describe the flutter point. Computations have also been done with the flow density fixed and the Mach number 
varied, but these results will be omitted for the sake of brevity. 

A. Linear Flat Plate Model 

7. Flutter of the Nextel in various configurations 

Table 6 includes the flutter densities of Nextel with in-plane loads, structural damping, and with a Pyrogel 
foundation. The results indicate that low-mode flutter is possible at M = 5 for a flow density significantly below that 
of air at sea level (p s i (20°C)=1.2 kg/m 3 ), for both simply supported and clamped boundaries on the Nextel alone. 
Addition of a relatively stiff Pyrogel foundation increases both the natural and flutter frequencies dramatically, 
which is likely non-physical due to the fact that the Nextel is not bonded at every point to the Pyrogel and the Kevlar 
motion is neglected. It is also seen that the addition of structural damping for the Nextel alone has little effect on the 
flutter results. This is consistent with the conclusions of Dowell. 3 Structural damping is much more significant when 
a Pyrogel spring foundation is added since the damping generalized force term is proportional to the fundamental 

natural frequency ( 4k! m ) . 


Table 6. Flutter results for the Nextel in various configurations, with M = 5. 


Configuration 

Flutter Frequency 
co f (rad/s) 

Flutter Density 

Pf(kg/m 3 ) 

Flutter 

Modes 

Nextel Alone, Nominal (BC: SS) 

2578 

0.51 

1,2 

Nextel Alone (BC: SS) 
Nx=Ny=1000 N/m 

2924 

0.61 

1,2 

Nextel Alone (BC: SS) 
Cn=0.05 

2606 

0.53 

1,2 

Nextel on Pyrogel Foundation 
(BC: SS) 

25679 

1.41 

1,2 

Nextel on Pyrogel Foundation 
(BC: SS) 

Cn=0.05 

27222 

5.72 

1,2 

Nextel Alone, Nominal (BC: Clamped) 

3954 

0.76 

1,2 

Nextel Alone (BC: Clamped) 
Nx=Ny=1000 N/m 

4224 

0.85 

1,2 

Nextel Alone (BC: Clamped) 
Cn=0.05 

3996 

0.8 

1,2 

Nextel on Pyrogel Foundation 
(BC: Clamped) 

25856 

1.67 

1,2 

Nextel on Pyrogel Foundation 
BC: Clamped 
Cn=0.05 

27443 

6.06 

1,2 
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Note that simply supported boundary conditions require that there be no deflection or bending moment at the 
edges. Clamped boundary conditions require that the plate slope must also vanish at the edges in addition to the 
deflection. Since this condition is difficult to enforce in an experimental setup, it is likely that the TPS test article is 
closer to simply supported than clamped. In any case, the Nextel alone with simply supported boundary conditions 
has the lowest flutter density of all the Nextel configurations. 


2. Flutter of the 3 -layer TPS 

The flutter results in Table 7 are for the 3 -layer TPS with various structural properties. Note that the boundary 
conditions indicated apply to the Nextel, Pyrogel, and Kevlar layers. 


Configuration of the 3- 
layer TPS 

Flutter Frequency eo f 
(rad/s) 

Flutter 
Density p f 
(kg/m 3 ) 

Flutter 

Modes 

Nominal (BC: SS) 

1326 

0.52 

1,2 

Pyrogel Mass = 0 
(BC: SS) 

2340 

0.52 

1,2 

N x =N y =1000 N/m 
(BC: SS) 

1662 

0.72 

1,2 

Cn= C,k = 0-05 
(BC: SS) 

1326 

0.52 

1,2 

Cn= C,k = 0.1 
BC: SS 

1339 

0.524 

1,2 

Nominal (BC: Clamped) 

2021 

0.76 

1,2 

Pyrogel Mass = 0 
(BC: Clamped) 

3579 

0.77 

1,2 

N x =N y =1000 N/m 
(BC: Clamped) 

2283 

0.93 

1,2 

Cn = C,k = 0.05 
(BC: Clamped) 

2026 

0.77 

1,2 

Cn= C,k = 0.1 
(BC: Clamped) 

2037 

0.79 

1,2 


Table 7: Flutter results for 3-layer TPS in various configurations. M=5. 


It is observed that low-mode flutter occurs for all configurations of the 3 -layer TPS. For any given configuration 
above, the clamped cases result in higher flutter frequencies and flutter densities than the equivalent simply 
supported cases. When the Pyrogel mass is neglected, the flutter frequency increases, which is expected since the 
effective inertia of the system is lower. Tensioning the TPS by a load of 1000 N/m only increases the flutter density 
by a factor of 1.2 to 1.4 depending on the boundary conditions. Since an in -plane load of 1000 N/m is quite large 
compared to what would be realistically applied, it is unlikely that tensioning will have a major effect on the flutter 
behavior in the wind-tunnel experiments. The effect of structural damping, even with a damping ratio as high as 0.1, 
also has little effect on the flutter densities. 

The dynamics of the problem may be further explored by comparing the results from the 3 -layer TPS to those of 
the Nextel alone. Consider the simply supported case only. The flutter densities are virtually the same for the 3 -layer 
TPS and Nextel alone, which further indicates that the Nextel stiffness is likely governing the flutter behavior. When 
the Pyrogel mass is neglected, the flutter frequency of the 3 -layer TPS approaches the flutter frequency of the Nextel 
alone, but is much lower when the Pyrogel mass is included since it adds inertia to the system. The effect of 
tensioning the 3-layer TPS by 1000 N/m has a greater effect than for the Nextel alone, which is likely due to the 
added stiffness of the bottom Kevlar layer that is introduced upon tensioning. Finally, the results in Table 7 illustrate 
the significance of allowing deflection of the Kevlar. The relatively large stiffness of the Pyrogel is included, yet the 
high frequencies and large flutter densities seen with the Nextel on Pyrogel foundation are no longer present. The 
freedom of the Kevlar now allows the Nextel to deflect without the Pyrogel restricting its motion. This is confirmed 
by relative motion seen in the flutter mode shapes computed with nonlinear theory in the next section. 
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Amplitude of Response 


B. Parabolic Curvature Model 

For this model, only simply- supported boundary conditions are considered, since the linear analysis indicated 
these conditions resulted in the lowest flutter density, and only the lower bound on flutter is desired. For the same 
reasons, structural damping and applied in-plane loads have been neglected in these calculations. Note that 
calculations with in-plane loading of 1000 N/m were performed, and the flutter density was unaffected at moderate 
to large H/h values. 


1. Flutter of the Nextel alone 

Before considering the effect of initial parabolic curvature on the flutter behavior of the Nextel alone, it may be 
of interest to compare the three dimensional nonlinear flutter results with H/h = 0 to the two-dimensional linear flat 
plate results computed earlier. Figure 1 1 shows that the peak in the flutter frequency spectrum of the nonlinear case 
is nearly aligned with linear flutter frequency. In addition, the flutter densities are the same, with p f ~0.5 at M=5. 
This consistency adds validity to the numerical methods and computational implementations for both models, 
though the possibility of a systematic error is not excluded. 

An example of the Nextel response just above the flutter boundary is shown in Figure 12 for H/h = 5. These limit 
cycle oscillations (LCOs) occur due to the geometric nonlinearity of the plate structure, which prevents the response 
from growing to infinity. It is evident that the response is highly chaotic, with average deflection amplitude of about 
five times the plate thickness. This result is consistent with the analysis of Dowell 4 , who found the flutter amplitude 
of curved plates to be on the order of the rise height, rather than the thickness. Also note that the baseline about 
which oscillation takes place is shifted below the zero point on the y-axis, which may indicate that plate initially 
deforms under the aerodynamic load and then flutters about that deformed position. 



© (rad/s) 

Figure 11. Flutter frequency spectrum for the Nextel 
alone with H/h = 0. The dotted line indicates the 
two-dimensional linear flutter frequency. 


4 


2 - 
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Figure 12. Nextel response (LCO) just above 
flutter boundary, with H/h=5. Time history is 
sampled at a fixed location on the plate. 


Table 8. Flutter densities for the Nextel alone at different values of H/h, with M = 5. 


H/h 

Curvature in X & Y 

Curvature in X only 

Flutter density p f (kg/m 3 ) 

# of streamwise modes 

Flutter density p f (kg/m 3 ) 

# of streamwise 
modes 

0 

0.5 

6 

0.5 

6 

1 

0.8 

6 

0.4 

6 

2 

1.1 

6 

0.4 

6 

3 

1.6 

8 

0.4 

8 

4 

1.6 

10 

0.4 

10 

5 

1.8 

12 

0.4 

12 
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The flutter densities of the Nextel alone as a function of H/h are given in Table 8. With streamwise and span wise 
curvature (X & Y), the flutter density increases as H/h increases. When curvature is present only in the streamwise 
direction (X), the flutter density generally decreases. A similar qualitative result for streamwise curvature was found 
by Dowell 4 with significantly different in-plane boundary conditions. The static aerodynamic loading due to the 
plate curvature (see Eqs. (1.16) and (A7)) is partially responsible for these effects, along with aerodynamic damping. 
Convergence is non-uniform and additional modes are required as H/h increases. The number of streamwise modes 
needed for convergence at each value of H/h is also given in Table 8. For these calculations, one span wise mode was 
retained. It was determined that in the case of only streamwise curvature, additional span wise modes have no effect 
on the flutter density. With spanwise curvature, the flutter density increases slightly with additional odd spanwise 
modes, but remains unaffected for additional even spanwise modes. Since the effect is minimal and usually within 
the iteration tolerance for the flutter density, only one spanwise mode was retained. We would expect spanwise 
modes to be critical if the flow were three dimensional, but the effect is likely less pronounced here due to the 
simplified aerodynamic theory used in this analysis. 

The flutter mode shape for the flat geometry is given in Figure 13, and the flutter mode shape for the 
parabolically curved geometry with H/h = 5 is given in Figure 14. The flat geometry exhibits low mode flutter 
behavior, which was also seen in the linear results. The amplitude near the flutter boundary for this case is on the 
order of the thickness. The behavior of the Nextel with parabolic curvature is chaotic, and there is some influence 
from higher modes. Since the chaotic motion doesn’t occur at a single frequency, there is no single flutter mode. The 
amplitude near the flutter boundary is on the order of the rise height, rather than the thickness. Analysis of the 
behavior in the full time series indicated snap -through buckling of the plate at non-periodic intervals. This behavior 
would likely be restricted in the wind tunnel layup configuration due to the presence of the rigid base. 



Figure 13. Flutter mode shape of the Nextel alone at an arbitrary 
time, with H/h = 0. Deflection has been normalized by the Nextel 
thickness. 



Figure 14. Flutter mode shape of the Nextel alone at an arbitrary time, 
with H/h = 5. The dotted line indicates the initial parabolic shape. 
Deflection has been normalized by the Nextel thickness. 
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2. Flutter of the 3 -layer TPS 

The flutter results for the 3 -layer TPS are given in Table 9. The flutter density at H/h = 0 is slightly larger than 
the density computed in the linear flat plate analysis, but by less than six percent. The flutter densities at all other 
values of H/h are slightly less than those of the Nextel alone. The number of streamwise modes needed for 
convergence at large H/h is also less than that of the Nextel. The limit cycle behavior is also chaotic for this system, 
similar to the behavior observed in Figure 12. Nextel response (LCO) just above flutter boundary, with H/h=5. Time 
history is sampled at a fixed location on the plate. Figure 12, but not shown here for brevity. 


Table 9. Flutter densities for the 3-layer TPS at different values of H/h, with M = 5. 


H/h 

Curvature in X & Y 

Curvature in X only 

Flutter density p f (kg/m 3 ) 

# of streamwise modes 

Flutter density p f 
(kg/m 3 ) 

# of streamwise 
modes 

0 

0.5 

4 

0.5 

4 

1 

0.8 

6 

0.4 

4 

2 

1.0 

6 

0.4 

6 

3 

1.2 

8 

0.3 

6 

4 

1.2 

10 

0.3 

8 

5 

1.4 

10 

0.3 

10 


The flutter mode shape for H/h = 0 is given in Figure 15. Similar to the Nextel alone, the flat geometry of the 3- 
layer TPS exhibits low mode flutter behavior, which was also indicated earlier by linear analysis. The high stiffness 
of the Pyrogel is evident because there is very little change in its nominal thickness during deflection. The free 
Kevlar base layer and Nextel generally deflect in unison along the panel length. Note that if the Pyrogel stiffness is 
decreased to well below its nominal value, it can be compressed easily and the Nextel and Kevlar are no longer 
strongly coupled. 

Unlike the flat geometry, chaotic, large amplitude flutter is observed for the parabolic geometry in Figure 16. 
Snap-through behavior was also observed in the time series for this case. 



x/a 


Figure 15. Flutter mode shape of the 3-layer TPS at an 
arbitrary time, with H/h = 0. Deflection has been 
normalized by the Nextel thickness. 


18 

American Institute of Aeronautics and Astronautics 



x/a 


Figure 16. Flutter mode shape of the 3-layer TPS at an arbitrary 
time, with H/h = 5. Deflection has been normalized by the Nextel 
thickness. 

3. Initial Conditions 

The simulations for both the Nextel alone and 3 -layer TPS do not require specification of initial conditions when 
H/h is greater than zero, due to the static pressure term in the aerodynamic forcing that provides a small initial 
deformation. However, these models respond quite differently when additional initial conditions are imposed. For 
example, when an initial displacement on the order of the thickness is imposed on the 3-layer TPS with H/h = 5, the 
flutter density drops by over fifty percent and the limit cycle behavior is initially less chaotic. Since the wind-tunnel 
test article is initially undeformed beyond the initial curvature, no additional initial conditions were specified for the 
results previously presented. 


C. Imposed Pressure Model 

As in the previous model, only simply- supported, stress-free boundary conditions are considered. While there is 
edge curvature present in both the streamwise and spanwise directions, it was determined that resolution of the 
spanwise shape was unnecessary for the flutter calculations, though this may not be the case for other types of 
calculations. As a result, only the first spanwise mode was retained. Flutter was computed for both the Nextel alone 
and the 3 -layer TPS. 

1 . Flutter of the Nextel alone 

The flutter densities for the Nextel alone are given in Table 10. The flutter density increases rapidly as H/h 
increases, which illustrates the inhibiting effect of static pressure on flutter. This effect was also observed by 
Dowell 8 for the case of a constant, uniformly distributed pressure on a plate. The number of streamwise modes 
required for convergence is also much larger than the parabolic curvature model at large H/h. This is because the 
additional stress from the imposed pressure results in coalescence flutter of higher natural modes, thus requiring 
these higher modes for convergence. 

Table 10: Flutter densities for the Nextel alone at different values of H/h, with M = 5. 


H/h 

Flutter density p f 
(kg/m 3 ) 

# of streamwise modes 

Peak in FFT 
(rad/s) 

Approx. Flutter 
Modes 

0 

0.5 

6 

2174 

1,2 

1 

1.9 

6 

6474 

3,4 

2 

4.0 

8 

10827 

4,5 

3 

6.6 

12 

20933 

6,7 

4 

9.0 

14 

26166 

7,8 

5 

11.0 

18 

48099 

9,10 
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Limit cycle oscillations (LCOs) of the Nextel just above the flutter boundary are shown in Figure 19 for H/h = 5. 
The period between t = Os and t = 0.005s is the time required to impose the static pressure, after which the 
aerodynamic pressure is applied. Despite a chaotic initial transient, the LCO is generally periodic, with amplitude on 
the order of the Nextel thickness rather than the rise height. This is in contrast to the LCO of the Nextel with 
parabolic curvature in Figure 12, which is chaotic with large amplitude oscillations near the flutter boundary. The 
periodic behavior of this system allowed for a peak frequency to be extracted from the Fourier transform of the time 
series results. These peak frequencies, assumed to be the flutter frequencies, are also given in Table 10. The 
approximate flutter modes were then determined by correlating the FFT peaks to the Nextel natural frequencies from 
linear theory in Table 3. 



Figure 19. Nextel response (LCO) just above the flutter 
boundary, with H/h = 5. The time history is sampled at a 
fixed location along the plate. Deflection has been 
normalized by the Nextel thickness. 

The flutter mode shape in Figure 20 illustrates two important features of the imposed pressure model. First, the 
leading edge of the Nextel is compressed from the oncoming flow, but the pressure-imposed geometry is generally 
retained. Just beyond the flutter boundary, the trailing edge begins to flutter at a relatively small amplitude, which is 
indicated by the red circle. 



x/a 


Figure 20. Flutter mode shape of the Nextel alone at an arbitrary time, 
with H/h = 5. The dotted line indicates the pressure-imposed shape. 
Deflection has been normalized by the Nextel thickness. 
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2. Flutter of the 3 -layer TPS. 

The flutter densities of the 3 -layer TPS are lower than the Nextel alone, which was also observed in the parabolic 
curvature model. These results along with the peak response frequencies are given in Table 11. The limit cycle 
oscillations are also periodic and similar to the behavior seen in Figure 19. While the flutter modes in this case are 
generally the same as the Nextel alone, the peak response frequencies are significantly lower. This is consistent with 
the lower natural frequencies for the 3 -layer TPS computed in the linear analysis. The flutter behavior in Figure 21 
is characterized by significant static deflection of the initial pressure-imposed geometry, especially at the leading 
edge, while the dynamic oscillatory behavior is generally found towards the trailing edge. 

Table 11: Flutter densities for the 3-layer TPS at different values of H/h, with M = 5. 


H/h 

Flutter density p f 

(kg/m 3 ) 

# of streamwise modes 

Peak in FFT 
(rad/s) 

Approx. Flutter 
Modes* 

0 

0.5 

4 

1450 

1,2 

1 

1.4 

6 

2695 

3,4 

2 

3.2 

8 

5413 

4,5 

3 

5.2 

12 

10466 

6,7 

4 

7.8 

14 

13035 

7,8 

5 

10.5 

16 

14959 

8,9 


*10 modes given to both the Nextel and Kevlar in the eigenvalue calculation 
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Figure 21. Flutter mode shape of the 3 -layer TPS at an arbitrary time, 
with H/h = 5. The dotted line indicates the pressure-imposed shape. 
Deflection has been normalized by the Nextel thickness. 


3. Initial Conditions 

As in the parabolic curvature model, no initial conditions are necessary in these simulations due to the initially 
imposed pressure. When additional initial conditions are specified, the computed flutter results are unchanged. It 
was also found that applying the static pressure and aerodynamic pressure at the same time (t = 0) produced the 
same results as delaying the aerodynamic pressure until the desired curvature was attained, as in Figure 19. 


VI. Conclusion 

The flutter behavior of a NASA HIAD thermal protection system in its wind-tunnel test configuration has been 
evaluated in this paper. Specifically, the proposed TPS configuration consisting of Nextel 440-BF20, Pyrogel 6650, 
and Kapton Kevlar Laminate was studied. Three theoretical structural models were developed based on 
photographic observations of the experimental test setup as described in previous conference proceedings 1,2 . The 
piston theory hypersonic flow approximation was used for the aerodynamic pressure in all cases. 
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The first structural model was a linear system of two-dimensional flat plates, which was solved using an 
eigenvalue analysis. With the Mach number fixed to a value of five, the flutter density parameters and flutter 
frequencies for various configurations were computed. It was determined that the Nextel outer layer governs the 
flutter behavior when coupled with the Pyrogel and Kevlar. The linear analysis also indicated that the TPS will 
flutter at Mach 5 with a flow density below that of air at sea level. Large in -plane loading of 1000 N/m applied at all 
edges slightly increases the flutter density parameter, however, more realistic loading would have a negligible effect. 
Structural damping was also found to have little effect on the flutter. 

The second structural model included the effect of initial parabolic curvature, which may arise from the restraint 
of the TPS test article at the boundaries. Von Karman’s large deflection theory was used, resulting in a nonlinear 
system of equations. The limit cycle oscillations of the TPS in this configuration are chaotic, with large flutter 
amplitudes on the order of the curvature rise height. It was determined that the flutter density could either increase 
or decrease with increasing rise height depending on the nature of the curvature. 

In the third structural model, curvature is imposed due to a static pressure from the rigid substructure, and the 
flutter behavior is significantly different. The effect of the static pressure is to increase significantly both the flutter 
frequency and flutter density, with periodic limit cycle response on the order of the Nextel thickness rather than the 
rise height. Plots of mode shapes indicated higher mode flutter at the trailing edge, rather than chaotic behavior seen 
in the parabolic curvature case. For both the parabolic and pressure-imposed curvature models, the 3-layer TPS has a 
noticeably lower flutter boundary than the Nextel alone for all nonzero rise heights. This indicates that the manner in 
which the TPS layers are stitched (or not stitched) will have an effect on the flutter behavior, but only if there is 
curvature present. 

Further examination of the TPS article in its experimental test configuration is needed to determine which 
structural model is the most realistic. In any case, the preceding analysis may help to suggest the best method for 
flutter suppression in the wind-tunnel tests. Out of all the models, the effect of imposing a static pressure on the 
underside of the TPS had the greatest effect in raising the flutter boundary, while in -plane loading and structural 
damping had less of an effect. The linear analysis also indicated that clamped boundary conditions had higher flutter 
densities than simply supported conditions. This suggests that modifying the restraint bracket in the test setup to 
prevent rotation at the TPS edges would also raise the flutter boundary. 


Appendix 

A. Modal Equations of Motion - Linear Flat Plate Model 

The modal equations resulting from the application of Lagrange’s equations are given in this section. All 
integrals have been written in dimensionless form for convenient numerical integration. The dimensionless variables 
are % = xl a and 7j = y / b . 
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The generalized forces, including aerodynamic pressure and structural damping are: 

Q n = ab[la n ^ EV £ £ p 2 dddp 
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Thus, the mass, damping, and stiffness matrices in the eigenvalue problem (1.10) are: 
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B. Modal Equations of Motion - Parabolic Curvature Model 

For brevity, the equations presented here are for a single (general) plate only, but the formulation can be easily 
extended to a 3 -layer system by following the approach used for the linear flat plate model (though an eigenvalue 
analysis is not possible here). Using Galerkin integration, the expression for the modal coordinate of the Airy stress 
function is: 
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The plate rise function has been non-dimensionalized using z = Hz . For convenience, let the large factor 
with Eh in the numerator be the quantity C(k) . Let the sum of the two integrals multiplying the product a n a t be 

I(n,i,k) and the sum of the three integrals multiplying a n be J(n,k ) . Thus, f k can be written as 


f k = C(k ) 
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This expression is then used in the expansion of the Airy function, along with the particular solution in Eq. 
(1.17). The generalized force due to aerodynamics and structural damping is: 
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The highlighted term in the generalized force expression is the static pressure loading due curvature, noted 
earlier as a critical factor in the flutter behavior of curved plates. 


The equation of motion for this system is: 
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The additional variables R(k ) , L(k) , and S(n) resulting from integral expressions in the stretching and bending 
energies are: 
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C. Modal Equations of Motion - Imposed Pressure Model 

The equations of motion for this case are also given by Eq. (A8), with the rise height H set to zero. This renders 
the plate initially flat, but the nonlinearity of the material is retained. Curvature will be imposed via a generalized 
force that corresponds to a pressure satisfying the static problem in Eq. (1.21). The expression for this generalized 
force is: 
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The mode shape functions a n are the same functions given by the expansion in Eq. (1.20) . Note that a single 
subscript is written in the equations of motion, but this implies summation over both streamwise and spanwise 
directions. The Fourier coefficients, expanded out in both directions, are: 


25 

American Institute of Aeronautics and Astronautics 



^ nm 


P\ 


(n n p x -2) cos (nnp x ) 2 p x cos (nnp x ) + 2 


3 3 

n n 


nn 


3 3 

n n 


H 


H \cos(nn p x )~ cos (nn p 2 )] 


nn 


+ 


+ 


H 


(n 2 n 2 -2)cos(nn) 2p 2 cos(nn) 2(p 2 +l) (2- n L n p 2 )cos(nnp 2 ) 2p 2 cos(nnp 2 ) 


2 2 2n 


3 3 

n n 


2(1 -2p 2 ) cos {nn p 2 ) 


■ + 


3 3 

n n 


Px 


nn 

2 2 2 


3 3 

n n 


+ - 


3 3 

n n 


- + - 


nn 


(n m p x -2)co s(mnp x ) 2p x cos(mnp x ) + 2 


3 3 

n m 


mn 


3 3 

n m 


+ — [cos(m^ p x )~ cos (mn p 2 )] + 
mn 


a-p 2 y i 


(nm - 2) cos (mn) 2 p 2 cos (mn) + 2 (p 2 + 1) 


3 3 

n m 


mn 


3 3 

n m 


+ (2-n m p )cos( W ^ 2 ) + 2p,cos(m^ 2 ) + 2(l z 2p ; )cos(m^ ; ) 


3 3 

n m 


mn 


3 3 

n m 


(A12) 


References 


! J. A. Del Corso, F. M. Cheatwood, W. E. Bruce, S. J. Hughes and A. M. Calomino, "Advanced High- 
Temperature Flexible TPS for Inflatable Aerodynamic Decelerators," 21st AIAA Aerodynamic Decelerator Systems 
Technology Conference and Seminar , Vol. 1, AIAA, Dublin, 2011, pp. 139-161. 

2 S. J. Hughes, J. S. Ware, J. A. Del Corso and R. A. Lugo, "Deployable Aeroshell Flexible Thermal Protection 
System Testing," 20th AIAA Aerodynamic Decelerator Systems Technology Conference , Vol l.,AIAA, Seattle, 
2009, pp. 291-311. 

3 E. H. Dowell, Aeroelasticity of Plates and Shells , Noordhoff International Publishing, Leyden, 1975, Chaps. 
2,3. 


4 E. H. Dowell, "Nonlinear Flutter of Curved Plates," AIAA Journal, vol. 7, no. 3, pp. 424-431, 1969. 

5 3M, "Ceramic Textiles and Composites Bulletin," 3M, Saint Paul, MN, 1996. 

6 The Dupont Company, "Kevlar Aramid Fiber Technical Guide," 30 October 2001. [Online]. Available: 
www.dupont.com. [Accessed 20 June 2012]. 

7 Microstructured Materials Group, "Silica Aerogels," Lawrence Berkeley National Laboratory, 1 January 2012. 
[Online]. Available: http://energy.lbl.gov/ecs/aerogels/sa-physical.html. [Accessed 20 June 2012]. 

8 E. H. Dowell, "Nonlinear Oscillations of a Fluttering Plate," AIAA Journal, vol. 4, no. 7, pp. 1267-1275, 1966. 


26 

American Institute of Aeronautics and Astronautics 



